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Abstract 

Starting from an Operator Product Expansion in the Heavy Quark 
Effective Tlieory up to order 1/ml we calculate the inclusive semilep- 
tonic decays of unpolarized bottom hadrons including lepton mass 
effects. We calculate the differential decay spectra dT/{dE-j-), and the 
total decay rate for B meson decays to final states containing a r 
lepton. 
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1 Introduction 



Recently, the theoretical investigation of inclusive semileptonic bottom hadron 
(Hb) decays has witnessed a renaissance through the use of the Operator 
Product Expansion (OPE) in the Heavy Quark Effective Theory (HQET) 
[H, ^. Up to now these investigations have been limited to the discussion of 
the zero lepton mass case. While this is a good approximation for inclusive 
semileptonic Hh decays involving an electron or a muon, a study of inclusive 
semileptonic Hh decays involving a r must certainly include the discussion of 
lepton mass effects. 

There are always two aspects to the lepton mass problem in semileptonic 
decays. First the decay dynamics changes as one is probing also the spin 
zero components of the exchanged gauge boson in addition to the spin one 
components active in the zero lepton mass case. In the formalism of semilep- 
tonic inclusive decays this means that there are two more structure functions 
(W4, W5) contributing to the decay as compared to the zero lepton mass case 
(1^1,1^2,^/3). In addition the decay kinematics gets changed and becomes 
somewhat more complicated than in the zero lepton mass case. As the very 
dynamics of the OPE method entangles the decay kinematics of particle and 
parton decay in a subtle manner one needs to carefully discuss the decay 
kinematics of the massive lepton case. 

Apart from theoretical considerations the calculation of inclusive semilep- 
tonic bottom decays has immediate phenomenological implications since cur- 
rent high energy experiments have reached a stage of sophistication that allow 
them to start measuring inclusive semileptonic Hi, rates involving the r. For 
example, recently the ALEPH collaboration reported on the first measure- 
ment of the inclusive semileptonic bottom hadron decay Hi, — > Xr^Ur [§]• 
They quote a branching ratio 

BriHb^XTUr) = (2.76 ±0.47 ±0.43)%. (1) 

In this presentation we shall closely follow the notation of , albeit with 
the necessary extensions to include lepton mass effects. We also include a 
separate section on the operator product method in HQET where we present 
general formulae which allow one to go beyond the 1/mf results considered 
in 0,§. 
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2 Operator Product Expansion in HQET 



Inclusive semi-leptonic B-decays are related to the transition operator 



where q{x) stands for up or charm quarks. The propagator term in 
be written in momentum space as 



1 — -1 
g(x)g(0) 



iS{x, 0) = i 
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The up or charm quarks may be considered as propagating in the external 
gluon field Fap generated by the heavy quark source. The expansion of the 
quark propagator in powers of the external field is best derived in the Fock- 
Schwinger gauge following the general discussion in 0]. One obtains 
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where Da = da + igA'^t"'. 

The OPE can only be meaningfully applied to (0) if short distance (large 
momentum) contributions are absorbed in the Wilson coefficients and the 
long distance (small momentum) contributions are absorbed in the matrix 
elements of the operators. We must explicitly extract the large mass scale 
into an exponential, i.e. pass to the operators in HQET. We write 



b{x) 



B(x). 



The field B{x) is related to the usual field by{x) of the HQET by 
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where D'^ = D^ — v^v ■ D. Using the expansion 



B{x) = B{0) + {x-D)B{0) + -{x-DfB{(}) + ... 
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where A;^ = iD^ relates to the residual momentum of the heavy quark, we 
obtain from (|[) 

V = /^/rf^xe-"(™^+'-^) -5(0)1^5(01,5(0). (8) 

The integrations in (^ are trivially carried out with the result 

V = BiO)T^S{mtV + k-q)T,B{0). (9) 

To obtain explicit results from this master formula we must substitute for the 
propagator in the background gluon field from (^) and expand the resulting 
expressions in powers of k using the formula 



1 
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p — m p — m p — m 

where p = m^v — q. 



... + ..., (10) 



3 Inclusive Decays Hi) X^^u + r + 

Consider a heavy hadron Hh moving with a four- velocity v and having a mass 
Mh^ (containing a b quark of mass mb) which decays weakly into a hadronic 
system of invariant mass Mx (containing a c or u quark of mass mj (j=c,u)) 
and an off-shell W~ with momentum q. The W~ decays subsequently into a 
r-lepton and an antineutrino z/^. Previous authors have investigated inclusive 
semileptonic Hi, decays involving the e- and /x-lepton where the mass of the 
final lepton can be neglected We extend the previous work to the case 

of the tau lepton where the r-mass (m,- =1.777 GeV) cannot be neglected. 
We denote the energies of the antineutrino and tau in the rest-frame of the 
initial hadron by E^, and Et-, respectively. 

The phase space in the variables (g^, Et-) has the form shown in Fig.l for 
the inclusive decay B XcTUj. and in Fig. 2 for the decay B — > X^TUr- -E-r 
may take values from mr up to (M^^ — {Mx^i„ + mr)^) / {2MhJ and is 
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bounded by 



< q < 



M2 + m2 - 2M^,,E, 



Here Mx^^^ = mo for b —>■ c and = m-Tr for b ^ u. 

The hadronic tensor of the inclusive process under consideration is defined 
to be 



X 



where Jj" is the hadronic current and j = u, c. The expansion of W'^'^ in 
terms of Lorentz invariant structure functions reads 

W^"" = -g'^'Wi + v^v^'W^ - ie^''<'''Vpq,W^ + q^q^'W^ + {(fv'' + v^q'')W^iYi) 

where q = Pr+Pu and Wi, .., are the structure functions to be determined 
from the dynamics. Note that there is no T-odd term proportional to {q^v'^ — 
v'^q'^) in the fully inclusive decay. 

The contraction of the hadronic with the leptonic tensor gives the spin 
averaged differential decay rate 



dq^dErdE^ 27t^ 

+Ws{q\E^ - E,) - ml{E, + E^)) 

+W, i^{q^ - ml)^ + W,{2mlE,)^ . (14) 

In comparison to the massless lepton case the factors which multiply Wi 
to change, and there are additional contributions from the structure 
functions and W^. 

The structure functions Wi, are given by the absorptive part of the 
matrix element of the transition operator 

T'^^ = {Hb\t^,\Hb) (15) 
= -g^'^'Ti + v^v''T2 - le^'^P'^ypq^T^ + q^'q^'T^ + {q^v"" + v^'q'')T^, 
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where ImT'^'^ = —nW^'^. What remains to be done is to determine the 
amphtudes Ti, ..,T^. 

First we expand the master formula of our operator product expansion 
in powers oi k. As we are interested only in operators up to dimension five 
we can terminate the expansion at order k"^. From the expansion of 



M7^(l - 75)(mfcy^- 4+ l^ + mj)Y{l-lb)u 
4((mf,f — g + kY — fn^j + ^e) ' 

in powers of k we obtain the contribution of the k^- terms to T^^: 

j.{0) ^ m - v.q 
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The T^\t^\t^'^ terms were already given in 0] . We list them here only for 
completeness. The contributions of the k^- terms to T^^ read: 

rrW T? ( ^ {mb-v.qf\ ^rrib , n \ ( ~^ , <f-i'"-<lf\ 

rpii) ( 1 2mb{mb-v.q) \ , '^^b . , ^ . ( ^ . 2m^\ 
= "^'"""[-Ao J+^^^^ + ^^Ha^ 

^ ^(^....)(|.) 

rr.il) p (mb-v.q\ 2mb ( -^h - v .q\ 

T5 = mbEb I — ^2 ) + ^(^^ + '^^) ( J ' ^ ^ 

where the spin-energy, kinetic energy and total energy matrix elements Gb, Kb 
and Eb, respectively, of the b-quark in the hadron Hb are defined as in M. 
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The order fc^ contributions are 

.(2) 



= ^ 
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There is an additional gluonic contribution originating from the gluon oper- 
ator in the second term of (H), i.e. 



^^6G"^e„^A.(mfe - qfig^'^Y + g^^^i^ - g^'^ + ^e'^'^^^7r]^^& (20) 
which leads to 



rp{g) _ _rr^Gb 

^3 - QA2 



3A^ 



Ti^) = 



(^) 



Summing the various contributions Tj = J2zTi , i G {0, 1,2,(7} we obtain 
, mbEb mlGb, . 
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T4 = f^iKb + Gb) 



where 



2Ao' 3 ' " 3A2 



-2{mb-v.q)mbEb 8ml Kb , 2 , \2\ , 2mlKb 



Ao 3A^ V. V . ^ 







Our formulae for Ti,T2,T3 coincide with the ones given in |Q, whereas the 
T4 and T5 expressions are new. 

We get the resulting differential decay rate ([T^ ) by taking the imaginary 
part of T'^'^ and dividing it by (— vr). This amounts to the replacements 



1 

a; 

1 

A^ 



5{{mbV - qf - mj), 





1 1 



Ag 



6'{{mbV - qf - mfj, 

5"{{mbV-qf -m]). (22) 



From the dispersion relation in the variable v.q it can be seen that this 
replacement is only legitimate if one integrates over the whole neutrino phase 
space. This is no restriction as the r-neutrino is not detected. 

As we are only interested in the double differential rate with respect to q^ 
and Er we integrate over Ey. Taking into account the (5-function structure 
according to (22) the integral over E^ takes the following form 

/""^ '""^^ dE4/i5(E, - El) + h5'{E, - El) + h5"{E, - E^)] , (23) 
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where 

The hmits of integration in (23) are as follows: i) for the values of {q^, Et) 
lying inside of region A of Fig.l (for which the minimum mass of the final 
hadronic state coincides with Mx^^^, the mass of the lightest hadron state 
containing a quark j), 

™m ^ q^ -ml 

2{E, + ^E^^ - ml) 

^•^ - ^r, (2bj 



2Mh, 



b 



ii) for the values of (g , Er) inside of region B of Fig.l, 



™ = (27) 

^rnax ^ Q - ^28) 

The continuous line which separates the two regions A and B is given by 
^ Ml^+ml-2MH,E, 

(29) 

As in the massless lepton case, the actual particle masses are so that 
always E^ < E^°^ as long as (g^, E^) stay within region A of the Dahtz plot. 
This is however no longer true for region B, where this inequality has to be 
imposed explicitly. Its effect is to give an upper bound on the mass of the 
final hadronic state (smaller than the kinematically allowed one (M^^^^ — m,-)). 
This is a new aspect which was not present in the massless lepton case. 

The net result of all the restrictions imposed by the 5 functions and its 
derivatives is that the double differential decay rate (23) is nonvanishing 
only within the partonic region delimited by the interrupted line in Fig.l, 
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determined in terms of quark masses. Its expression is 

dV 
Tiydydq^ 

e(^+)e(-^_) (l2(-g^ + 2q^y + q'p - fr] - f - y^ - yp + yrj + y + pr) - r}) 

+12E,{2q^ - 2fy - 2fp + 2fr) + yp-yri + y- 2pr]) 

+SKh{-q'^ + fp - fr] + 2f - 3y + pr] + 677) 

+8(76(2?^ - 2fy - 2fp + 2fri - f + yp-yr)-2y- 2pr] + 5r))) 

-5{z+) [2AE,a+i-q^y - 2gV+ + 2f + y^ + 2ya+ - yrj - 2y - 2a+r] + 2r]) 

+4:Kb{fy^ + 4:fya+ + 4gV+ - 20gV+ - Afr] - y^ + y^r] + Aya\ + Ay a +7] 

+12ya+ + Ayr] + Aa^r] - 20(7+77 - 
+16G6(^^ - fy(T+ - fy - 2^V+ + fa+ -f + y^cr+ + 2yal - ya+r] + ya+ 

+yr] — 2a1r] - a^r] - rf + 77)) 
+5{z-) {2AEb(T^{-q^y - 2gV_ + 2q^ + + 2ya- - yr] - 2y - 2a^r] + 2r]) 

+AKb{fy^ + Afya_ + 4gV^ - 20gV_ - 4g^7; - y^ + y^V + Ayat + Aya_ri 

+ 12ya_ + Ayr] + Aa'^rj - 20(T_7/ - Ar]"^) 
+16Gb{q'^ — (f'yo'- — (f'y — 2(f'a^_ + (fa_ — + Z/^ct- + 2ya^ — ya^r] + ya_ 
+yr] — 2(7^7/ — (T_77 — 77^ + 7^)) 

+%5\z+)a+Kb (Aq^ - qY - - - + 4g'?7 + Z/^ + 4z/V+ 

-y^V + %cr+ - 4?/(7+77 - 4(7+77) 

-8S'{z_)a_Kb (4g^ - fy'^ - Afya_ - Afy - Afa'i + Afr] + y^ + Ay'^a_ 

-y'^r] + Aya'i - Aya-r] - Aa'ir]^ (30) 

where 

z± = l + f-p-y- 2a±, (7± = {f - r])i±, ^± = / = , (31) 

y ± Vy^ - 477 

rrib ml m^ m^ 
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and 



^ 1927r3 ■ 

In order to obtain the lepton spectrum, the double differential decay rate 
(30) has to be integrated over q^. The integration limits are 



^1 



2{l + r]-y) 



(33) 



The upper limit corresponds to 2;+ = and the lower one (f_ to z_ = 0. 
The result is 

dr 



Tbdy 



y2 _ 



2(— + 3yri + 3y — Ai]) — Gyp 



6{-4:T] + y + yr]) 
{1 + rj-yy 



P 



-16?7 + 6y + Grjy — 2y'^ 



4(-y2 + 3y- 2r]) 



-P 



12(— Gr] + 21]"^ + y + 3i]y — rjy"^ 



{l + r]-yy 
4(20?7 - 4:1]^ -6y - I4:r]y + 4y^ + 4:r]y^ - y^) 



(i + v-yY 

_l(_20r/ + 9y + 2y2) 



4(-20?7 + 32r7^ + Ar]-^ + 3y + 2T]y - 25r]^y + 2y^ + Ar]y^ + 2r]^y^ - 2y^ + rjy 

{l + ri-yY 

^ (80r/ - IQAif + 8r/^ - 18y - mrjy + bOrfy + lOy^ 



3(l + r/-y)5 



--(-20r7 + 15y + 2y^ 



8(4?7 - 3y - 3r7?/ + 2?/^) 4(12?7 + 477^ - 3|/ - 15?7y + 3y^ + 2?7?/2) 2 



(1 +7]-?/)^ 



(1 +r]-y)'^ 



8(-20r7 + 4772 + 6?/ + 14?7y - 4y2 - 4r^y2 + y ^ ^ 
3(l + r7-y)4 ^ 



(34) 
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The total decay rate can be obtained by integrating the spectrum formula 
over the partonic range of values for y 

2^<y<l-p + V, (35) 

with the result 

r = r,[/o(p, v) + EJe{p, v) + k^p, v) + GJg{p, v)] ■ (36) 

Here 

fo{p,v) = VR{l-7ri-7ri' + ri^-7p+12rip-7ri^p-7p^-7rip^ + p') 

-24^1 - /) log ^ + + (p - v)] (37) 

fE{p,ri) = /R(5 - 1977 + 577^ - 377^ - 19p - 477P + 2I77V + + 21?7p2 - 3p3-) 

-24[v\l + 3/) log ^^^"^""^ + (p - v)] (38) 
/i^(p,77) = 2v^(-3 + 1377 + 77^ + 77^ + 13p - 477P - 777V + - 777p^ + p^) 

+48[77^(1 + p2) log ^ + ^-p-^ + ^ rj)] (39) 

/g(p, 77) = 2\/R{-l + 7r] + 7rf-rf + 7p-l2r]p + 7rfp + 7p^ + 7r]p^-p^) 

+48[r7^(l - p2) log ^ + + (p - 77)] (40) 



where 



i? - [1 + 7]- pf - At]. (41) 



The limit of a massless final state quark p — requires a special treatment 
because of the appearance of 5{y — 1 — 77) and 5'{y — 1 — rj) functions in the 
lepton spectrum. They are due to the fact that in this limit z+ vanishes 
along the entire line y — 1 + r]. By considering as a function of y at fixed 
g^, one can write 
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Inserting these relations into Eq.(30) and performing the g^-integration be- 
tween the hmits 

ql^\{y±^f^^) , (44) 

one obtains 

— - = 7y2_4^ r2(_2y2 + + 3^ _ 4^) 
+4E,(-y2 + 3y-2n+^{l- r] fd{y - 1 - v)) 
-^K,{-20rj + 9y + 2y' + (1 - r))\5 - 2ri)S(y - 1 - v)) 

~G,{-20ri + 15y + 2y^ - (1 - 77)^(4 - ri)S{y - 1 - v)) 
+8K,{l-r))h{y,ri)5'{y-l-rj). (45) 
Here h{y, rf) is defined by 



Miiy) (1 - r)i-i +^+y + 2^+(r -^)] 

(4g^ - g^y^ - 4:fya+ - Afy - Afal + Afr] + y^ (46) 
+4y V+ - y^77 + Ayal - Aya+r] - 4(7+7?) 

and satisfies 

hm —h{y,v) = --{l-7i)\ (47) 

An integration of the spectrum (45) over y within the hmits 2y^ < y < 1 + rj 
gives the total rate 

r = r.lgoiv) + E,gE{v) + K,gK{v) + Gbgoir])] (48) 

with 

gQ(r]) = 1 - 87] + 8r]^ - T]"^ - 12r]'^ logr/ (49) 

gj^{r]) = 5 -2Ar] + 2Ar]'^ -8r]^ + 3r]* -Uri'^logT] (50) 

gj^{ri) = -6 + 32r/-24r/2-2r/^ + 24r/Mogr/ (51) 

g^^r]) = -2 + 167] - 167]^ + 2r;^ + 24r/^ log 7] . (52) 
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This is the same as what one would obtain by directly taking the limit p — 
in (36). In doing the ^/-integration of (45), the end-point y = 1 + rj must be 
considered as included in the region of integration. 

Before plotting the lepton energy spectrum we need to specify the values 
of the nonperturbative matrix elements Gb and Kb. The matrix element Eh is 
given by Ef, = Gb + Kb- is related to the hyperfine splitting of degenerate 
heavy meson partners and is taken as Gb = —0.0065 as in [Q. For Kb we 
take Kb = 0.01. Further parameters are rub = 4.8 GeV and rric = 1.39 GeV. 
We shall plot the energy spectrum for mesonic b-decays into c and u. 

Let us first discuss the decay spectrum (34) for b — > c plotted in Fig.3. 
The leading order contribution corresponding to free quark decay is larger 
than the 0(1/ ml) corrected spectrum over most of the y-range. The correc- 
tions become larger as the lepton energy increases. The corrected spectrum 
takes an upward turn to the high end of the spectrum. The OPE can, how- 
ever, no longer be trusted in local sense at the high end of the spectrum. As 
a rough estimate where this happens we demand that the nonperturbative 
terms should not exceed 30% of the leading parton model result. Judging 
from the curve in Fig.3 this means that we trust the energy spectrum lo- 
cally up to 1/^0.98. However, assuming duality, the integrated spectrum is 
expected to be given correctly even including the high end of the spectrum. 

The total inclusive decay rate B XcTUr is 

r/Ffc = (0.122 + 1.286 ■ Eb - 1.409 ■ Kb - 0.245 ■ Gb) = 0.122 - 0.008 (53) 

for rrir = 1.78 GeV, rrij = rric = 1.39 GeV and rrib = 4.8 GeV, whereas for 
rrir = we would have had 

F^^=o/r6 = (0.543 + 3.361 ■ Eb - 3.904 ■ Kb - 1.086 ■ Gb) = 0.543 - 0.020 . 

(54) 

The lepton mass effects appear thus to increase the procentual contribution 
of the nonperturbative I /ml corrections in the total result from -3.7% to 
about -6.6%. 

Next turn to the tau spectrum for b ^ u, which is plotted in Fig.4. Also 
shown is the free quark decay spectrum for comparison. Again one should 
trust the spectrum locally only for y^O.98. 

The total decay rate is 

T/Tb = (0.371 + 2.584 ■ Eb - 2.954 ■ Kb - 0.741 ■ Gb) = 0.371 - 0.016 , (55) 
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where the value niu = has been taken. 

While preparing this manuscript for publication we became aware of the 
paper of Koyrakh |TD[ , who has also considered lepton mass effects within the 
operator product expansion method. While we agree with his final numerical 
results, we do not agree with all his formulas, which we suspect to be marred 
with misprints. 
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Figure Captions 



Fig.l The Dalitz plot for the inclusive decay B — > XcTUr- The interrupted 
line shows the region of the phase space which is populated in the 
parton model. 

Fig. 2 The Dalitz plot for the inclusive decay B — > XuTUr- The interrupted 
line shows the region of the phase space which is populated in the 
parton model. 

Fig. 3 T-energy spectrum for inclusive scmilcptonic bottom meson decays 
6 — c, normalized through division by Ti,. The continuous line is 
the spectrum including 0(l/m^) corrections and the dotted line shows 
the leading parton model result. 

Fig.4 T-energy spectrum for inclusive semileptonic bottom meson decays 
6 — > M, normalized through division by F;,. The continuous line is 
the spectrum including 0(1 /ml) corrections and the dotted hne shows 
the leading parton model result. 
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